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Abstract 

We consider the method of self-similar renormalization for calculating critical 
temperatures and critical indices. A new optimized variant of the method for 
an effective summation of asymptotic series is suggested and illustrated by 
several different examples. The advantage of the method is in combining 
simplicity with high accuracy. 
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I. INTRODUCTION 



In the theory of critical phenomena one usually obtains critical temperatures and critical 
indices as expansions in powers of some parameters which, as a rule, are not small. For 
instance, in the Wilson e— expansion [1] one has e = 1. In the field-theoretical approach [2] 
the expansion is in powers of the renormalized coupling constant g = 1.4. Such expansions, 
as is well known, are asymptotic and lead to reasonable results only in the low orders. The 
direct use of higher orders makes the results only worse [2,3]. When a number of terms 
in a divergent series is known one may invoke resummation techniques, such as Pade and 
Pade-Borel ones. While the knowledge of only a few first terms does not permit to use 
these techniques. Thus, one always confronts the problem of how to improve the results of 
a divergent series having only a few terms. 

In the present paper we suggest a solution of this problem advancing a method which 
has the following peculiarities: (i) It permits to accomplish an effective summation of a 
divergent series consisting of just a couple of terms, when no other method is applicable; 
(ii) It is simple and accurate, providing an accuracy not worse than sophisticated Pade 
and Pade-Borel techniques involving about ten terms, when they are available; (iii) It is 
regular, unambiguously prescribing the way of action. The method suggested is a variant 
of the method of self-similar renormalization [4-6]. The latter is a renormalization-group 
approach using self-similarity of subsequent perturbative terms. Since renormalization group 
is nothing but a kind of a dynamical system, the approach can be formulated in the language 
of dynamical theory with the usage of its powerful techniques as well as of those of control 
theory [6-8]. Then the number of the approximation order plays the role of discrete time. 
Motion with respect to the latter corresponds to transfer from one approximation to another. 
This makes it possible to define a dynamical system whose trajectory is bijective to the 
sequence of approximations. Such a dynamical system with discrete time has been called 
the approximation cascade, Convergence of a sequence of approximations is equivalent to 
stability of a dynamical trajectory. The stability and, respectively, convergence are governed 
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by control functions. The fixed point of a trajectory defines the sought function. 

This method has been successfully applied to the eigenvalue problem in quantum me- 
chanics [9,10]. Here we consider the resummation problem for the asymptotic series in the 
theory of critical phenomena and advance a novel variant of the method, more appropriate 
for this problem. 



The complete description of the method with the corresponding mathematical foundation 
can be found in Refs. [4-8] . In this Section we give only the general scheme of the method, 
which is necessary for explaining the new variant we suggest. 

Suppose we are interested in a function f(x) of the variable x G (—00,00). Let this 
function satisfy a complicated equation that cannot be solved exactly. Assume that by 
means of perturbation theory we can get a sequence {pk{%)} of perturbative approximations 
Pk{x), where k = 0,1,2..., enumerates the approximation order. Usually, perturbation se- 
quences are divergent. To extract a meaningful result from a divergent sequence one has to 
involve the so-called resummation techniques. In the method of self-similar renormalization 
a divergent sequence can be made convergent by introducing additional functions governing 
convergence (see [4-8]). These functions, because of their role are called governing or con- 
trol functions. Let s be a set of such control functions entering into a sequence {Fk(x, s)} 
obtained by a perturbation algorithm. 

In addition to introducing the control functions, the main idea of the method of self- 
similar renormalization is to treat the passage from one approximation to another as a 
motion with respect to the approximation number k — 0, 1,2,... This motion is realized in 
the functional space of the considered function as follows. Define the initial approximation 



as an equation for the expansion function x = x(f, s). Substitute the latter back to Fk, so 



II. METHOD OF SELF-SIMILAR RENORMALIZATION 




(1) 



that 
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y k (f,s) = F k (x(f,s),s). (2) 

The relation inverse to (2) is 

F k (x, s) = y k (F (x, s), s). (3) 

Let {y k } form a group of transformations with respect to k — 0, 1, 2... Then the trajec- 
tory {y k (f,s)} of this dynamical system, according to definitions (2) and (3), is bijective, 
that is, in one-to-one correspondence, to the approximation sequence {F k (x, s)}. This dy- 
namical system with discrete time k has been called [7,8] the approximation cascade. The 
attracting fixed point of the cascade trajectory is, by construction, bijective to the limit of 
the approximation sequence {F k (x, s)}, that is, corresponds to the sought function. 

To deal with continuous time is easier than with discrete. Therefore, we embed the 
approximation cascade {y k } into an approximation flow {y(t,...)} with continuous time 
t > 0. This implies that the trajectory {y(t, f, s)} of the flow passes when t — k — 0, 1, 2, 
through all the points of the cascade trajectory, 

y(k,f,s) = y k (f,s) (A; = 0,1,2,...). (4) 

The evolution equation for the flow reads 

^y(t,f,s) = v(y(t,f,s)), (5) 

with the right-hand side being the velocity field. The latter, in the language of 
renormalization-group theory, is often called the Gell-Mann-Low or (3— function. 

Integrating the evolution equation (5) from t = k to t = k* , we get the evolution integral 

I**, dj_ = k . _ kt 

Jy k v(f, s) 

in which y k = y{k, /, s) and y k+ i = y{k*, /, s). Before specifying the numbers k and k* in 
the limits of the evolution integral, let us note that the differential form (5) of the evolution 
equation, or its integral form (6), are equivalent to the functional relation 

y(t + t',f,s) = y(t,y(t'J, S ),s). (7) 



The latter in physical applications is labeled as the self-similarity relation, which explains the 
term we use. The self-similarity, in general, can occur with respect to motion over different 
parameters. In our case, this is the motion over the steps of a calculational procedure, the 
number of steps playing the role of effective time. 

If there exists an attractive fixed point of the approximation-flow trajectory, then it is 
always possible to find a number k* in the evolution integral (6) such that the upper limit 
yl would correspond to an expression 

F*(x,s) = y(k*,F (x,s),s) (8) 

representing, with the desired accuracy, the sought function f(x). If yl would be an exact 
fixed point, then (8) would give an exact answer to the problem. However, a fixed point 
can be reached only after infinite number of steps k — > oo. For a finite number k, the limit 
yl may represent a fixed point approximately, because of which it is named the quasi -fixed 
point. Our aim is to reach the latter as fast as possible, that is, during the minimal time 

t* k — min(A;* - k), (9) 

or the minimal number of steps. When there are no additional restrictions, the minimal 
number of steps counted by k is 1, so that 

abs min t* k — 1. (10) 

In the case when some constraints are imposed on the motion, the minimal time (9) should 
correspond to the conditional minimum. For instance, if a value / = f(x ) of the sought 
function f(x) is given for some £ , then we can find t* k by requiring the trajectory of the 
approximation cascade to pass through the given point / . . 

To calculate the evolution integral (6), we need to define the velocity field. This can be 
done by the Euler discretization of (5) yielding the finite-difference form 

Vk(f,s)=y k (f,s)-y k - 1 (f,s). (11) 
Substituting (11) into (6), and using (3), we come to the representation 
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for the evolution integral (6), where F k = Fk(x, s), F£ +1 = F£ +1 (x, s). 

Finally, we have to define the set s of control functions. The role of the latter is to govern 
the convergence of the approximation sequence. This convergence can be expressed, in the 
language of dynamical theory, as the stability of the cascade trajectory. A useful tool for 
analyzing stability is a set {[x k } of the local multipliers 

d 

Vk(f,s) = —y k (f,s). (13) 

The inequality 

I /»*(/,*) I<1 (14) 

is the condition of local stability at the step k with respect to the variation of an initial point 
/. The equality | /i k (f,s) \— 1 implies local neutral stability. For a convergent sequence 
corresponding to a contracting mapping, the condition of asymptotic stability is 

IM/>*)hO (fc->oo). (15) 

The approximation cascade {y k } describes the motion in the functional space {/}. To 
return to the domain {x}, we must use the inverse transformation (3). With the help of the 
latter we may pass from the multiplier (13) given on the functional space {/} to its image 

m k (x,s) = ii k (F (x,s),s) (16) 

being a function of x. For the image (16), the same stability condition as in (14) can be 
written, 

| m k (x, s) \< I. (17) 

According to (15), the local multipliers diminish when approaching an attracting fixed point. 
That is, the variation of initial condition / produces lesser and lesser effect on the trajectory 
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as soon as the attractor becomes closer. In other words, the lesser are the absolute values 
of multipliers, the more stable is the trajectory. Therefore, it is reasonable to define the 
control functions as those minimizing the absolute values of the local multipliers, making by 
this the trajectory more stable at each step k. In this way, a set s of control functions is to 
be defined by the principle of maximal stability written as 

| m k (x, s k (x)) |= min s | m k (x, s) \ . (18) 

Because of this, the control functions s k (x) defined by the principle (18) may be called the 
stabilizing functions or stabilizers. 

Note that the control functions may be introduced in several ways, as is discussed in 
Refs.[4-10], however always being related to stability conditions and the closeness of a tra- 
jectory to an attracting fixed point. In all the cases the control functions are to be defined 
so that they could accomplish their main job, i.e., to govern the convergence of an approx- 
imation sequence, which, in the terms of dynamical theory is equivalent to stabilizing the 
cascade trajectory. In the present paper we shall use the definition of stabilizers given in 
(18). 

After the stabilizers are defined, we have to substitute them into the corresponding 
approximations F k (x, s) getting 

f k (x) = F k {x,s k (x)). (19) 

This stage can be called the stabilizing renormalization of a perturbative sequence. 

Then, considering the motion near the renormalized quantity (19) by means of the evo- 
lution integral (12), we obtain 

fi(x) = FZ(x,8 k (x)). (20) 

This step can be called the dynamical renormalization. And the whole procedure of the 
double renormalization (19) and (20) is named the self-similar renormalization. It is worth 
noting that the evolution equation (5) is, generally, nonlinear and can have several different 
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solutions leading to different self-similar approximations (20). In such a case, to select a 
physically meaningful solution, we need to involve additional conditions as constraints. The 
role of the latter can be played, e.g., by properties of symmetry, asymptotic properties at 
x — > or x — > oo, sum rules or other relations containing some known information on the 
character of the sought solution. Such additional constraints narrow the set of possible 
solutions to a class with desired properties. Thus, we should always remember from what 
class we are looking for a solution. 

Keeping in mind that we wish to get a good accuracy for the sought function, having just 
a few perturbative terms available, we need to find out some tricks which could effectively 
increase perturbation order. We suggest below one such trick. 

Suppose that there is a sequence of approximations Pk{x) having polynomial structure, 
k showing the order of the polynomial. This order can be effectively increased by means of 
the multiplicative transformation 

P k {x,s) =x s Pk {x), s > 0. (21) 

Then, the order of the expression (21) becomes k + s. The transformation inverse to (21), 
as is evident, is 

p k (x) =x~ s P k (x,s). (22) 

Following the method described above, we consider the sequence {P k (x,s)} and construct 
an approximation cascade whose trajectory {yk(f, s)} is bijective to {Pk(x, s)}. Solving 
the evolution integral (12), we have P£(x,s). From the principle of maximal stability (18) 
we define the stabilizers Sk(x). Substituting these into P£(x,s) and invoking the inverse 
transformation (22), we obtain the self-similar approximation 

ft(x)=x->M PZ{x,8 k (*))- (23) 

The multiplicative transformation (21) looks as the most natural for the case when the 
perturbative approximations Pk(x) have the form of polynomials or series, generally speaking, 
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not necessarily in integer powers. The factor x s effectively increases the approximation order, 
and s plays simultaneously the role of stabilizer. 

What power s we have to choose, that is, to what effective order we need to go is 
dictated by the principle of maximal stability selecting the most stable trajectory of the 
approximations cascade. In particular, it may happen that s — 0, and we do not need to 
proceed further, or, vice versa, we may have to go to the limit of s — > oo, thus making 
allowance for all approximation orders. In each concrete case, an effective order which we 
need to reach depends on how good is the perturbative sequence {pk(x)} we start with and, 
respectively, how much information can be extracted from its first terms by means of the 
double renormalization (19) and (20). 

The optimization by introducing the stabilizing control functions into the powers of 
perturbative polynomials renormalizes the algebraic structure of the latter. Because of 
this, and in order to distinguish the suggested optimization procedure from other possible 
variants, we shall call it the algebraic self-similar renormalization. 

To concretize the procedure, let us write explicitly 



k 




(24) 



n=0 



as a polynomial of the order k. Following (21) define 



k 



Pk(x,s) = J2 



,n-\-s 



(25) 



n=0 



Similarly to (1), we have 



P (x,s) =a x s = f, 



(26) 



from where the expansion function is 




(27) 



The definition (2) yields the points 




n=0 



(28) 
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of the approximation-cascade trajectory. For the velocity field (11) we get 

^+i(/,s) = a fe+ i(-)^ +1 . (29) 
ao 

From the evolution integral (12) we find 

P k+1 = ~ 777^ — Z k+i s ■ ( 30 ) 



I (fc+1) a k+1 tt p — ^ 

1 1+r-: — r k ) 



k + l 



„ -+ 1 

s a n 



The multiplier (13) becomes 



M/,s) = E^(i + ^)(-f)^ (si) 

n=0 °0 S a 



and its image (16) reads 



m k (x,s) = J2-a + 1 )x n - (32) 

The principle of maximal stability (18) defines the stabilizers s k (x), whose explicit expres- 
sions depend on the coefficients a n . According to the transformations (21)-(23), from (30) 
we obtain 



Pk{x) 



Jk+i — ffc+ii a,,, *? , , ,k+i,_s_-> y^) 



where Sk(x) defines the most stable trajectory. When there are no additional conditions, 
the minimal value t* k — 1, as in (10). 

As is noted above, it may happen that the most stable trajectory corresponds to s — > oo. 
Let us show how the self-similar approximation (33) simplifies in this case. It is straightfor- 
ward to check that the limit of the right side in (33), as s — > oo, leads to 

r k+1 (x)=p k (x)eM—x k+1 )- (34) 

a 

One may notice that renormalizing p k (x) in (34) we can obtain the recurrence relation 

f; +1 (x) = r k (x)eM^x k+1 )- (35) 
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It is possible also to derive several other relations permitting to repeat the self-similar 
renormalization several times, which is useful when working with high-order terms. However, 
in what follows we shall limit the consideration of particular examples by keeping only a 
few terms of the corresponding perturbative series. This is to emphasize that the method 
suggested allows to reach good accuracy with a minimal number of perturbative terms, when 
no other resummation technique is applicable. Comparing (33) with (34), we see that the 
self-similar renormalization can yield quite different expressions, from the fractional form to 
exponential one. Below we shall illustrate this by some simple examples. Each appearing 
form of an approximation results from choosing the most stable trajectory by which it is 
possible to reach a quasi-fixed point during the minimal time. Recall in this connection the 
analogy with classical mechanics. Notice also that it is possible to follow a trajectory which 
is stabilized by imposing additional conditions, such as asymptotic properties, or prescribing 
that the trajectory is to pass through some given points. In such a case the motion will not, 
generally, be accomplished during the absolute minimal time (10), but the latter should be 
defined from the additional constraints imposed. All these variants will be exemplified in 
the following sections. 



with < x C 1. How could one continue this expression from small i < 1 to x > 1, when 
no other information is available? 

Following the algebraic self-similar renormalization, we construct the transformed poly- 
nomial (21) or (25) , which for the case (36) is 



III. ILLUSTRATION BY SIMPLE EXAMPLES 



Suppose that by perturbation theory we have got 



= 1 - x, 



(36) 



Pi(x, s) 



= x - X 



(37) 



According to (26), we have the expansion function 
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x(f,s) = fV s . (38) 

Then, Eq.(28) gives 

yi(f:s) = f-f 1+1 *. (39) 

The velocity field (29) becomes 

v 1 (f : s) = -f 1+l *. (40) 
The evolution integral (12), leading to (30) now yields 



p*( x ,s) = (^y. (4i) 



For the multiplier (32) we have 



m±(x,s) = 1 - (1 + -) x. (42) 
s 

Minimizing the absolute value of the latter gives the stabilizer 

s J x ) = 0<x<l, 

' l ~ x ~ ~ (43) 

si(x) — > oo, x > 1. 
The self-similar approximation (33) reduces to 

Being interested in the region x > 1 we have to take the limit s\(x) — > oo. Therefore, for 
the self-similar approximation (44) we obtain 

f*(x) = lim (— — ) s = e"*. (45) 

In the same way, a linear expansion 

= a + aix, (46) 

with a , ai 7^ 0, derived for | x |<^C 1, can be continued to the region | x \>\ ^ | where it is 
represented by the self-similar approximation 
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f*(x) = a exp(— x). (47) 

do 

Thus, we may conclude that the exponential (47) is a general self-similar representation of 
a linear expansion (46), when no additional constraints are imposed. 

Now turn to the case, when we want to construct a self-similar continuation of (36) 
satisfying the prescribed asymptotic behavior 

f(x) oc x n , x — > oo, (48) 

where n > is fixed. Repeating the same procedure as earlier, we come to (44). Comparing 
the latter with (48), we get Si(x) = n, so that 

n(x) = <^r. (49) 

Generalizing this result for a linear combination (46) under the asymptotic condition (48), 
we have 

/i*(*)=«o(— ^ ) n (50) 

na — aix 

In this way, one perturbative expansion may have several self-similar representations cor- 
responding to different imposed constraints. The form of these representations can vary 
between the exponential, (47), and fractional one, (50). However, for each given constraint 
this form is uniquely defined. If no constraints are imposed, the form of the resulting self- 
similar approximation is governed by the stabilizers obtained from the principle of maximal 
stability of a self-similar trajectory. 

Let us illustrate how accurate is a self-similar approximation and how it is possible 
to increase the accuracy by considering higher-order terms of a perturbative expansion. 
For this purpose take the function ln(l + x) with x > 0. We opt for this function, as 
an example, since the logarithmic expressions are typical of thermodynamic potentials in 
statistical mechanics and of generating functionals in field theory. 

Write down the three first perturbative approximations for ln(l + x) in powers of 
x, thinking that x is small: 
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pi(x) = X, 

p 2 (x)=x-^ ) (51) 
Ps(x) = x - 4 + t- 

Our aim is to construct self-similar approximations for ln(l + x) in the region x ~ 1, with 
expansions (51) that are valid only for x <C I. 

Following the standard prescription of the method, define the transformed polynomials 
(25) for those in (51), which gives 



Pi(x, s) = x 



1+s 

1 



P 2 (x, s )=x 1+S -^, (52) 
P 3 (x,s) = x 1+s -^ + ^. 

Now, the initial, i.e, the lowest order approximation is Pi(x,s). Therefore, as in (26), from 
the equation Pi(x, s) — f we find the expansion function 

x(f,s) = f&. (53) 

The points of the approximation -cascade trajectory (28) are 

yi(f, s) = f, 



!&(/,*) = vi(/,*) ( 54 ) 



2+£ 

2 

y 3 (f, s) = y 2 (f, s) + | /£f. 



3 

For the velocity field (21) we get 

v 2 (f,s) = -\ 
vs(f,s) = I 

The evolution-integral solutions in (30) become 



(55) 



P ;(x, s) = [ 2 } 1 + s ^ x ]i+s P *( x s) = [ 1 r^] 1 ^- (56) 

2(l + s)+x J 3V ; [ a; 2( 1 _|)Tfi 3(l + s) J v ; 

And for the multipliers in (32) we have 
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m 2 (x,s) = l-§(f±f) x, 

(57) 

m 3 (x, s) = m 2 (x, s) + |(fjf) x 2 . 

The stabilizers are to be defined at each step by minimizing the absolute values of the 
corresponding multipliers in (57). For instance 

s 2 ( x ) = 0, < x < 1, 
s 2 (x)=2£t, l<x<2, 

S2(x) — OO, X > 2. 

The corresponding expression for S3(x) is also easy to find. For the region of interest, 
where x ~ 1, we have 

s 2 (x) = 0, s 3 (x) > oo (x = 1). 
This leads to the self-similar approximations 

= (58) 
/•(x)=x(l-f)exp(^), 

obtained from (56) as in (23). 

In order to check the accuracy of (58) as compared to the perturbative expansions (51), 
define the percentage errors 

e k (x) ee Pk{ *l~ f l x) x 100% (59) 

and , respectively, 

= ^[J^ x 100% ' ( 6 °) 

where f(x) = ln(l + x). At the point x — 1, the errors (59) and (60), calculated with respect 
to /(l) = ln(2) = 0.693, are 

e 2 (l) = -28%, e 3 (l) = 20%, 
e*(l) = -3.8%, 4(1) = 0.67%. 
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As is seen, the accuracy of the self-similar approximations in (58) is an order higher than 
that of the perturbative expansions in (51), and this accuracy can be increased by taking 
into account additional perturbative terms. 

In this section we considered simple examples in order to make transparent all steps of 
our method. This will permit us in the following sections to avoid the repetition of the 
technical details when applying the method to more complicated physical problems. 

IV. CALCULATION OF CRITICAL TEMPERATURE FOR 2D AND 3D ISING 



In this Section we calculate the critical temperature T c of the two-dimensional (2d) and 
three-dimensional (3d) Ising model starting from approximate expressions for T c obtained 
by the variational-cumulant expansion (VCE)[11]. The convergence of VCE approximations 
is very slow and even using the seven consecutive approximations one obtain T c for 2d Ising 
model with the percentage error of about 11%. We use below the simple variant of the self- 
similar renormalization when the exact value of the sought function is known for some point, 
namely the knowledge of T c for the 2d Ising model will be used as an optimization condition 
for the trajectory, determining the optimal number of steps t*, then used for calculating T c 
for the 3c? Ising model. Although the expressions below are a little complicated for getting 
the result in an explicit form, we can realize here the numerical variant of the self-similar 
renormalization , when sought function is obtained implicitly. 

We rewrite the expressions for the critical temperature from [11] in terms of a new 
variable x — \ , where d is the space dimension. Then we calculate T c (x), related to T c as 
follows: T c = \%(x). 

Write down the three approximations to T c (x), 



MODEL 



T cl (x) = 2-x, 



12-12x+2 
6-3x ' 



(61) 




24-36a:+83: 2 +5a: 3 
12-12S+2X 2 ' 
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which correspond to the second, third and fourth approximations of [11], respectively. 

Following the standard approach described above in Section 3, from the equation T c (x) = 
f we find the expansion function x(f) = 2 — f. The points of the approximation cascade 
trajectory are 

Viif) = f, 

VM = f (^), (62) 

For the velocity field we get 

v (f) = l (19/ 2 -22/ + 4)(/-2) 2 
3UJ 6 (/2 + 2 /-2) / • 

The evolution integral cannot be calculated explicitly, so we expressed T* 3 implicitly as 
a function of x and t* given by the equation 

r^a* df 



r^ 3 _dj_ = f 



'Tc2 V 3 (f) 

and obtained T c * 3 numerically for two different values t* = 1, or t* = 1.5515, corresponding 
to a non-optimized and optimized variant respectively. The latter, optimizing number of 
steps, was obtained from the condition T c * 3 (|,t*) = where T 2 d = 2.269 is the celebrated 
Onsager temperature for the 2d Ising model. 

At t* — 1 and d = 2 we found that T* = 2.531 and the percentage error of our estimate 
is eg(T c ) = 11.542%, approaching the percentage error reached using seven consecutive 
approximations to T c [11]. 

We should point out here that the authors of [11] did not attempt to calculate T c for 
3c? Ising model, where the best known "exact" numerical value of the critical temperature 
is T c = 4.51 [12]. We obtained T* = 4.712 at t* = 1 and d = 3. The percentage error is 
eg(T c ) = 4.47%. For the optimized t* = 1.5515 our estimate of T* = 4.548 with the error 
equal to 0.838% is quite accurate. 

It is worth noting that, in principle, the parameter x = h may be considered as "small" 
for d = 3 and the expressions for T c (a;)can be expanded in powers of x, thus presenting T c 
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in the form of " i— expansion" . We performed the same renormalization procedure as above 
for 

Tcz{x) ~ 2 — x — ^| 

and found that in 2d for t* = 1, T c * = 2.453 with e*(T c ) = 8.109% and in 3d-case with 
t* = 1, we have T c * = 4.701 with e*(T c ) = 8.109%. 

For the optimized t* = 1.241 and d = 3, T c * = 4.624 and e^(T c ) = 2.527%. We see 
that the renormalized —expansion gives by order of magnitude the same accuracy that is 
reached from the renormalized original expressions. 

One can conclude from the results presented in this Section, that both the rate of con- 
vergence and accuracy of the VCE are greatly improved by applying the self-similar renor- 
malization to the starting VCE approximations. Situation encountered while renormalizing 
T c was somewhat very simple, since we possessed three reasonable approximations for renor- 
malization and also knew the exact value of T c for the 2d Ising model. In the next Section 
we meet the case when the number of the terms available are not sufficient for any variant 
of renormalization discussed above and no exact value for the quantity under consideration 
is known. 

V. RENORMALIZATION OF SHIFTED-POWER EXPANSIONS 

In this Section we apply a modified variant of the self-similar renormalization to 
the " shifted-power expansions" for the critical indices [13]. Shifted-power expansion or 
A— expansion is of particular interest for systems without upper critical dimensionality. It 
can be applied to the calculation of critical exponents of a system described by a Landau- 
Ginzburg (LG) Hamiltonian. Within the framework of one of many possible realizations of 
the method, the leading nonlinear term of the LG Hamiltonian,-?/' 4 , is replaced by (-?/> 2 ) 3 ~ A , 
then A is used as an expansion parameter and at the end of calculations one should set 
A = 1. The expressions for critical indices 77 and v were obtained in the following form: 
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r] = b 2 {n) A 2 , 
v = a + a 2 (n) A 2 , 



(n+4)(re+2) 
48(3n+22) 2 ' 



(63) 



ao = \, a 2 (n) 



(n+4)(ra+2) 
12(3n+22) 2 ' 



where n is the number of the components of the order parameter. 

The corrections to the mean-field values given by (63) are about two orders of magnitude 
too small. The authors of [13] had noticed that the results are strongly influenced by the 
if> 6 interaction and remain too distant from the analyzed if) 4 behavior. This conclusion 
agrees well with the rigorous results of [14] where it was shown that special Gaussian points 
n = —4, —2 should appear when the if) 6 model is considered, while for the ip 4 model only the 
point n = — 2 exists. Appearance in the expressions (63) of the combination n + 4 signalizes 
that the critical indices in the interesting physical region n = 0, 1...3 are influenced by the 
point that does not have any meaning for the if> 4 model at all. Nevertheless, the Gaussian 
point n = — 2 does have physical meaning for the if) 4 model [15-17] and, formally, the 
results given by (63) are correct at n — —2. Therefore, one may hope that a systematic 
and pernicious influence of the point n = —4 can be weakened by some renormalization 
procedure, at least when a physical region not very distant from the point n = — 2 is 
considered. We should stress also that the expressions in (63) do not obey another rigorously 
studied limit of the LG Hamiltonian, when n — > oo, i.e. r)(n — > oo) ^ and u{n — > oo) ^ 1 
[18]. Therefore we do not expect that a successful renormalization of (63) could be realized 
for n very distant from the point n = —2. 

Our approach to the renormalization of r\ and v should vary, since for the index v the 
two terms in the A— expansion are available and the renormalization could be carried on 
straightforwardly, while for the index rj only single term was obtained and it is not possible 
to proceed with extra assumptions. 

Write down the two consecutive approximations to the index v in powers of A : 



and apply the same procedure that leads to the expression (47) with a substitution of x to 



v (n) = a , 
v 2 (n) = a + a 2 (n) A 2 , 
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A 2 . We also retain in the final expression for the renormalized critical index z/| an effective 
time t* which will be exploited as an optimization parameter: 

u*(A,n,t*) = a exp{^^A 2 t*}. (64) 

Setting here A = 1 and t* = 1 we see that no considerable effect was achieved and the index 
v remains in the interval v = 0.5 to 0.509 while n varies from n = —2 to infinity. Impose 
now an additional condition that in the case of n — 0, corresponding to the random walk 
problem [19], the approximation cascade trajectory pass through the value of the critical 
index v = 0.588 known approximately, but with very high accuracy for this physical problem 
[20]. The choice of this point for optimization of the trajectory is dictated also by the desire 
to receive renormalized values for the physically interesting region n — 1,2,3 of "true" phase 
transitions using information only from the physically distant region, where the random walk 
problem may be a good choice because it does not correspond to a "true" phase transition, 
but only a formal analogy exists with the n — > limit of the LG Hamiltonian. 

So, from the condition v%(l,0,t*) = 0.588 we obtain t* = 59 and the following values 
at n — 1, 2; 3 : 

1, = 0.633, 1/2(1,2,0 = 0.676, z/ 2 *(l,3,t*) = 0.715. 

These values are reasonable as compared to the experiment, high-temperature series, and 
Borel-summation results [ 21]. We should point out that only by means of a single parameter 
we obtain simultaneously renormalized values for the physically interesting situation, i.e 
a systematic multiplicative error in the initial expansion, caused by peculiarities of the 
A— expansion can be eliminated by a single renormalization step. 

It is also possible to find optimal t* from the condition restoring the correct value of 
v at n — > 00, where, as is shown rigorously, v — 1. In our case, the results happen to 
be much better than for the initial A— expansion, but still largely underestimate v, e.g at 
n — 0, v* — 0.544. Nevertheless , this variant of optimization is of general interest because 
often the limit n — > 00 (or d — > 00) is well known. 
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In the case of the index r] one should proceed differently, since only a single term in 
the A— expansion is known: we added to the expression for r\ the term linear in A with 
some yet unknown positive coefficient b defining thus a new quantity fj. Then carried out the 
renormalization procedure for fj, repeating the steps leading to Eq.35. From the renormalized 
quantity fj* using the variational condition = 0, we determine b as a function of n, A, t*. 
Finally, we subtracted the term b A from fj* to find rj*. Following this prescription define 

fj = bA + b 2 (n)A 2 , 

then find 

^* = 6Aexp{^4^A t*}. 
b 

This yields 

v * = fj* -bA = t*b 2 (n)(e - 1) A 2 . 

At t* — 1 the results still remain too small. Imposing an additional optimizing condition 
on t* by analogy with the case of v, that rj* should equal to 0.026, where the critical index 
rj for the random walk problem is taken from [20], we find that t* = 44 and 

rj*(n = 1, t*) = 0.038, r]*(n = 2, t*) = 0.048, r]*(n = 3, t*) = 0.057. 

These results overestimate r], especially for n — 3, but are much more realistic than the 
initial value rj ~ 10~ 4 — 10~ 3 . This systematic error can be understood if to notice that already 
the initial A— expansion does not obey the limit of n — > oo, and 77(00) = 0.002 instead of 
zero. This systematic deviation cannot be fully corrected by a variational-renormalization 
procedure. The same is true in the case of index u, but in this case more information is 
available from the initial A— expansion and the results of renormalization remain reasonable 
even at n — 3. 

We conclude that by applying the self-similar renormalization to the A— expansion for 
the critical indices, one can obtain reasonable estimates for rj and u, although further im- 
provement of these estimates does not seem plausible, since initial expressions violate an 
exact relation in the n — > 00 limit and possess an unphysical Gaussian point at n — —4. 
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In the next Section we meet the case when both limits at n = — 2 and n — > oo are 
violated. 



We have seen in the previous Section that the A— expansion mimicking the widely ac- 
cepted Wilson e— expansion is, in the best case, a crude approximation for the critical indices, 
since an important n — > oo limit is violated already in the starting terms of the A— expansion. 
The question naturally arises whether the e— expansion obeys exact limits for critical indices, 
namely at n = — 2 and n — > oo ? The discussion of this question for all critical indices 
will be presented later. In this Section, we consider the Wilson e— expansion for the critical 
index 5, discuss its n = — 2 and n — > oo limits, observe that they are violated and suggest 
the self-similar renormalization approach allowing to restore the correct limiting values for 
the index 5. 

Consider the well-known Wilson e— expansion for the critical index 5 [1,18] up to the 
quadratic terms in e : 



At n = — 2 and n — > oo, 8 — 4.5. From the exact results for n = — 2 vector model [16] 
and from the scaling law for the d = 3 the result 5 = 5 follows. The same value 5 = 5 was 
obtained in the case of the spherical model (n — > oo, d — 3) [18]. The percentage error for 
the critical index 5 in these limits is therefore —10%. In the physical region n — 0, 1...3 the 
e— expansion gives 5 = 4.47 — 4.46. Unfortunately, there is no much experimental data or 
theoretical estimates available for the index 5, but if we accept the scaling laws as correct 
and estimate from these values of 5 another critical index r] = |=| , which is much better 
known, both experimentally and theoretically, then we appear in the physical region with 
7] ~ 0.1, that largely (by three times) overestimates r\. One may think that the index 5 
is underestimated by the Wilson e— expansion. So, we need by means of the self-similar 



VI. IMPORTANCE OF ASYMPTOTIC PROPERTIES 




(65) 
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renormalization procedure to continue the asymptotic expression (65) to the region of e ~ 1 
with simultaneous restoration of the incorrect limiting values at n = — 2 and n — > oo. 
Introduce the following consecutive approximations to the quantity 5 = 5 — 3: 



Si(e) = e, 
62(e) = e + c 2 (n)e 2 

By repeating the same steps that led us to (56), we obtain 



6 U+s 



(66) 



5 2 (e,s) = [— 

and 



1+s 



5* = 3+[— e -^] 1+s . (67) 

1 l+s 

Setting in (67) e = 1 and n = — 2 , or n — > 00, it is easy to show that only for s = both 
limits can be satisfied! Therefore 

5* = 3 + - . (68) 

1 - c 2 {n)e 

The expression (68) in the physical region gives the following values: 

5*(n = 0) = 4.882, 5*(n = 1) = 4.862, 
5*(n = 2) = 4.852, 5*(n = 3) = 4.847. 

The index 77, corresponding to these values can be easily obtained from the scaling law: 

v ( n = 0) = 0.02, r](n = 1) = 0.024, 
r)( n = 2) = 0.025, r?(n = 3) = 0.026. 

These results better agree with the general understanding of 77 as of a small index and 
are much closer to the results of the Borel summation than the initial rj = 0.1, obtained 
from the e— expansion (65) and the scaling law. 

We conclude that the application of the self-similar renormalization improves the Wilson 
e— expansion for the index 5 both qualitatively and quantitatively. This example stresses 
once again an importance of obeying different reasonable limits in the expressions for the 
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critical indices. Another illustration is given in the next Section where the self-similar 
renormalization of ^—expansion and of the e— expansion around lower critical dimension 
two (in powers of d — 2) is considered. 



VII. INVERSE LARGE COMPONENT EXPANSION (1/N) AND EXPANSION 

IN POWERS OF D-2 (2+e). 

The large n— expansions (^—expansion) [18,22] and e— expansion around the lower crit- 
ical dimension two (d — 2— expansion) [23] had raised high expectations as an alternative 
to the Wilson e— expansion and field-theoretical approach [2]. Nevertheless they had never 
became competitive, remaining a useful guide to the region of large n and a good qualitative 
tool, when different aspects of the two-dimensional behavior are considered. It is clear that 
the values of critical indices given by - and d — 2— expansion do not obey the n = —2 
Gaussian limit, becoming divergent at n = and n = 2 (d — 3), respectively. Therefore, 
it is not possible to get a reasonable estimate for v and r\ in the physical region n = 1, 2, 3 
lying too close to the spurious pole and too far from the correct n — > oo limit, supported 
by both expansions. 

Consider the —expansion for the critical index 7 : 

24 

7 = 2-—, 

T\ A n 

from where 

7 (n = 1) = -0.432, 7(71 = 2) = 0.784, 7 (n = 3) = 1.189. 

Correspondingly, the two approximations in powers of \ can be written down 

7oM = 2, 
7l (n) =2--^. 

Proceeding in accordance with the self-similar renormalization prescriptions, the multi- 
plier mi(n, s) can be found: 
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As is seen, the minimum of | mi(n, s) | for n > 2 is reached for s — > oo. This gives 

7*(n) = 2exp( — \-), 

and, correspondingly, 

7 *( n = 2) = 1.089, 7l *(n = 3) = 1.334. 

We see that for n = 3 the value given by 7* becomes reasonable, deviating from the result 
of the Borel-summation 7 = 1.386 [20] with the percentage error equal to —3.752%, while 
the initial ^-—expansion has the percentage error of —14.214%. For n — 1, the minimum of 
I mi(n, s) I is reached for s = 5.633; correspondingly 7*(n = 1) = 0.508. We observe that 
for n = 1, 2 the results are improved if compared to the initial -^-—expansion . 

The ^—expansion for critical index rj is given as follows [18]: 

This equation becomes singular at n = and negative at n — —2, so that, despite its correct 
by design behavior at n — > 00, the values given by (70) at n — 1, 2, 3 are too large: 

77(71 = 1) = 0.076, ??(n = 2) = 0.086, r?(n = 3) = 0.068. 

The direct application of the self-similar renormalization using ^ as a renormalization 
parameter with 

7] - 8 1 7] - 8 1 ( 8 ) 3 1 1 
37r 2 n' 37r 2 n 3 7r 4 n 2 ' 

as consecutive approximations does not improve the situation, since the influence of singu- 
larity at n = is too strong. To avoid this divergence, we re-expanded the expression (70) 
in powers of the parameter y = fej^p , expressing n as a function of y : 

n = ^(l-16 2 / + V / l-24y). (71) 
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This choice of the re-expansion parameter is not unique, but the combination fre- 
quently appears in the Wilson e— expansion. Such a transformation restores the correct 
value of 77 at n = — 2 and also keeps intact the correct limit at n — > 00. Up to the third 
order in y we obtain: 

, o , 8 , 112 1856 14336 

7i = ay + by+cy, a = — , b = c = — - — . 

Thus the following consecutive approximations may be written down: 

m(y) = a v, 
m(y) = ay + by 2 , 

Vs(y) = ay + by 2 + cy 3 . 
Proceeding in the usual manner, we obtain 

nKv) = ayexp (-y), (72) 

a 

r ] ;(y) = (ay + by 2 )eM-y 2 ). (73) 



a 

Returning to the initial variable we have 

r/ 2 *( n = 0) = 0.013, r/ 2 *(n = 1) = 0.016, 
rj*(n = 2) = 0.018, r)* 2 (n = 3) = 0.019, 

and 

r)*( n = 0) = 0.015, r/*(n = 1) = 0.02, 
77*(n = 2) = 0.023, r/|(n = 3) = 0.025. 
The values of multipliers in these cases are as follows: 

m 3 (y,s) =mi(y,s) + ff±f y 2 , 

These values are very close to each other, e.g for n — 3, m 2 = 1.022 and m 3 = 1.037. 
From the stability viewpoint the corresponding approximations are almost equivalent. It is 
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also possible to improve results for 77 by applying the second self-similar renormalization, as 
in the recurrence relation (35). This gives 77* in the form of a continued exponential [24]: 

77* = ayexp(-yexp(^y)), 
a 

so that 

77* (n = 0) = 0.016, 77*(n = 1) = 0.024, 
rf(n = 2) = 0.03, 77*(n = 3) = 0.032. 
These values, especially for n — 2, 3, are quite reasonable. The results for n — 0, not 
surprisingly, remain too small, since we used for the renormalization procedure only the 
large n expansion, obviously too short of information about the limit for small n. In order 
to weaken the influence of the particular way of defining coefficients a, b, c, it is possible 
to proceed with a variational-optimization procedure, considering 77* as a function of two 
unknown parameters a, b and determining them from the conditions 

^=0, *£ = 0. 



For the particular choice 



we find that 77* = 4eq/ 3 and 



da db 



77* = ayexp(-y - _ £ ^ ), 



rf( n = 0) = 0.019, r}*(n = 1) = 0.032, 
v *(n = 2) = 0.04, 77*(n = 3) = 0.043. 

The e— expansion with e = d — 2 around the lower critical dimension for the critical index 
u, is written in the form 

z/ -i =c 7_ 2 + ^-^! (74) 

for d > 2 and n > 2. At n = 3, d = 3, 1/ = 0.5, giving rather crude estimate coinciding 
with the mean-field result. The self-similar renormalization using e as a parameter for 
renormalization does not improve this value, because v — > 00 in the starting point d = 
2, n — 3. On the other hand, at n — > 00, d = 3, formula (74) gives 1/ = 1, i.e the correct 
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limiting value known rigorously for the spherical model. Re-expanding (74) in powers of ^ 
around this correct value we obtain 



1 „ 1 



v- y = 1 + - + 2— + ... (75) 
n n 2 

Proceeding in accordance with the standard prescription and using ^ as a renormalization 
parameter, we define 

^ = 1, 
^\n) = l+ 1 -. 
Then, we readily obtain for v*(n) the following expression 

vl(n) = exp(--) 

and i>l(n = 3) = 0.717, giving reasonable estimate for the critical index v. The percentage 
error, when compared to the result of the Borel summation v = 0.705 [20] equals 1.702%, 
while for the initial d — 2— expansion it equals —29.078%. 

For the critical index rj, the d — 2— expansion has the following form [23]: 

1 n — 1 

v = a (n)(d - 2) - b(n)(d - 2) 2 \ a(n) = -, b(n) 



n-2' v ' (n-2) 2 ' 

At n = 3, d = 3, this gives r\ = —1, in disagreement with all known about this index. 

We use below d — 2 = e as a renormalization parameter, since at d — 2, r\ = 0, being a 
reasonable starting point for the trajectory. Define the following approximations to i] : 

?7i(e,n) = a(n)e, 
r) 2 (e, n) = ?7i(e, n) - b(n)e 2 . 

The multiplier m 2 (e,n) = 1 — ^)f-j-fe reaches its minimum at s — > oo, therefore 

r] 2 (e,n) = a(n)eexp{ --e), 

a[n) 

and 772(1, 3) = 0.135, which is a considerable improvement if compared to rj = — 1 from the 
d — 2-expansion. 
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We conclude that the self-similar renormalization improves the quality of estimates also 
for - and d — 2-expansions, achieving the quantitative agreement with other approaches. 
However, the broken n = — 2 Gaussian limit still makes the possibilities of improving the 
results very narrow, usually an improvement is achieved for n = 3, but not for lower n. 

In the next Section we consider an opposite case of an expansion obeying only n = — 2 
limit, but with broken limit at n — > oo. Such situation is similar to that of Section 5, but 
no non-physical Gaussian points will be present. 

VIII. SELF- AVOIDING WALK PROBLEM (N+2-EXPANSION FOR N=0) 

Interesting properties of the LG model for n = —2 have been analyzed in a number of 
works [14-17]. Physically, n = — 2 corresponds to a Gaussian polymer with the exponents 
7 = 1, rj = 0, v = 1/2. From the scaling laws in d — 3 one can see that a — \, 8 — 5, (3 — \. 
It seems natural to develop expansions in powers of n + 2 (n + 2— expansion) around this 
well defined limit [16]. To our knowledge, this idea has never been put into practice. The 
task of obtaining the n + 2— expansion is simplified if we note that the Wilson e— expansion 
for the critical indices 7, 77, u, a, f3 obeys the n = — 2 limiting values. In order to obtain the 
n+ 2— expansion we simply re-expand the Wilson e— expansion at e = 1 in powers of n + 2. Of 
course, n + 2— expansion could be derived also from the "first principles" in a way similar to 
the e— expansion , or —expansion. The nearest to the point n = — 2 physically interesting 
case is located at n — 0, corresponding to the self-avoiding walk problem equivalent to 
a polymer. We believe that the case of the order parameter dimensionality of n = — 2 
and n = are closely connected in a way similar to the connection existing between the 
space dimensionalities d — 4 and d — 3, with the only difference that the e— expansion is 
substituted by the n + 2— expansion. We apply below the self-similar renormalization to the 
n + 2— expansion for the critical indices, presenting only the results for n — 0. The values 
of the critical indices are not as good for n = 1,2,3. The expansion parameter is too large 
in the latter cases and also the n — > 00 limit is violated, so that the trajectory strongly 
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deviates for larger n from the correct but distant starting point. 

Up to the second order in e, the critical index v is given as follows: 



1 n + 2 e n + 2 2 2 ^ 

2 4n + 8 8n + «^ V 



Put here e = 1 and expand (76) in powers of z = n + 2, up to the third order terms in z : 

1 „2 , * 5 1 7 



f = o + + 02-2 + d^z ; ai = — , a 2 = —7377, 03 = 



2 J ' 96 864' 3456 

The following approximations to the quantity v = v — \ can be readily written down: 

V\{z) = d!Z, 

V2{z) = aiz + a 2 z 2 , (77) 
£3(2) = a x z + a 2 z 2 + a 3 z 3 . 

Following the standard way we find the multipliers 

m 2 (z,s) = 1 + ^^z, 

ai 1+s ' 

m 3 (z,s) =m 2 (z,s) + f 2 ^ s z 2 . 

Both multipliers at the point z — 2 {n — 0), reach their minimum at s = 0, where 
m 2 (2,0) = 0.911, m 3 (2,0) = 0.444. Consequently the trajectory restored using all three 
approximations from (77) will be more stable, than that restored from only two approxima- 
tions, both trajectories being stable. The evolution integral (12) gives 

v 2 (z) = a lZ - — , 1/3 {z) = — 2 1 » 

l ~t z {l-2^v 2 {z))^ 

and v 2 (z = 2) = 0.6, v%{z = 2) = 0.588. The former value is exactly the Flory "mean- 
field" exponent [25], and the latter is the same as the Borel-summation result, considered 
as the best known estimate for polymers [20]. The latter value z/| = 0.588 should be 
trusted more, since it is obtained moving along the more stable trajectory than v\. It 
is encouraging that the Flory and field-theory results, in our consideration come out as 
successive approximations. 

The Wilson e— expansion for the critical index (3, up to the second order in e, is 
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H 2 2(n + 8) 2(n + 8) 3 V ; 

So, the n + 2— expansion, up to the third order, becomes 

15 1 1 



f3 = a + a x z + a 2 z + a 3 z , 



4' A 144' ' 864' ° 432 
This results in the following approximations for (3 : 

Po(z) = a , 

Pl(z) = Po(z) +CHZ, 

(80) 

(3 2 (z) = ^(z) + a 2 z 2 , 
(h(z) = [3 2 (z) + a 3 z 3 . 

The multipliers m,i(z,s), m 2 (z,s), m 3 (z,s) for z = 2 reach their minima at s — > oo; so 
that 777,3(2, 00) < m 2 (2, 00) < mi (2, 00). The evolution-integral solution for (3 3 (z) is 

P;(z)=p 2 (z)eM-z 3 ), (81) 
ao 

and = 2) = 0.301. This coincides with the result of the Borel summation [20]. 

The critical index 77 has the following e— expansion, up to the third order: 

and the corresponding n + 2— expansion, up to the third order terms, can be obtained : 

2 , 25 23 43 

77 = a\z + a 2 z + a 3 z ; ai = — — , a 2 = — , a 3 



864' * 2592' 31104 
Thus, the following approximations result 

7/1(2) = aiz, 
7] 2 (z) = a\z + a 2 z 2 , 
7/3(2:) = a x z + a 2 z 2 + a 3 z 3 . 
The multipliers m 2 (z,s) and 7773(2;, s) at 2: = 2 satisfy the condition | m 2 (z, s) |< 1, | 
7773(2;, s) |< 1 for arbitrary s, the former satisfying the condition | 7772(2;, s) \— at s — 0.586, 
the latter becoming minimal at s = 0. Thus, we find 

*/ N 1 */ n 



i-^-) 1+s MW (1-2^(^)1 
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with 772(2 = 2) = 0.034, 773(2; = 2) = 0.023. The former value is very close to the so- 
called unconstrained e— expansion (77 = 0.031 ± 3 [26]) and the constrained e— expansion 
(77 = 0.0320 ± 25 [26]). The latter value approaches closely the result of Borel summation, 
77 = 0.027 ± 4 [26]. The scaling law \v{l + 77) = j3 is ideally satisfied with our v\ = 
0.588, r]* = 0.023 and (3* = 0.301. 

It is also worth noting that a non-optimal, but stable trajectory for 773(2,5 — > 00) = 
772(2) exp(^-z) leads us to the point 7/3(2, 00) = 0.027, which is exactly the value of the Borel 
summation. 

The e— expansion for the index 7, up to the second order, has the form 

n+2 n+2 /9 n9 . . 

~> = 1 + ^8f + 4t^W in+22n + ° 2)t - (83) 

and the n + 2— expansion, up to the third order, is 

7 



7 = 1 + a x z + a 3 z 3 ; a x = — , a 3 



72' J 216 
The multiplier at z = 2 acquires its minimum value at s = and 

^( Z ) = Q _ 2o3 z 2W2 + X 

with 73(2 = 2) = 1.165, and with the percentage error 63(7) = 0.345% as compared to the 
results of Borel summation [26]. Moving along a non-optimal but stable trajectory with 
s — > 00, we come to 73(2,5 — > 00) = a!2exp(^-2 2 ) + 1, which yields the value 73(2,3 — > 
00) = 1.161, in the complete agreement with the Borel summation [20,26]. The scaling law 
i/(2 — 77) = 7 gives for 1/ = i/3, 77 = 773 the value 7 = 1.162. The percentage error in this case 
equals 0.258%. 

For the critical index a, the e— expansion, up to the second order, reads: 

a = A - n e - ^ + 2)2 (n + 28)e 2 , (84) 
2(n + 8) 4(n + 8) sV ; ' K J 

and the n + 2— expansion, up to the third order, follows: 

1 .2 3 1 1 1 



a — — — a\z — a,2Z — azz , ai = -, a,2 = 7^7:, ^3 



2 6' 432' 108 
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For the quantity a = ^ — a, with the set of approximations 

ai(z) = CLiZ, 
a 2 {z) = ai(z) + a 2 z 2 , 
a z (z) = a 2 (z) + a 3 z 3 , 

following the conventional prescriptions, we find the solution corresponding to the most 
stable trajectory with s = : 

—*/ x oi{z) 

013 {Z) [l-2(a 3 /a?)cT 2 (z)]V2' 

and al(z = 2) = 0.217. The field-theory Borel summation results are not available. From 
the scaling law find that the value v = 0.594, corresponding to ctg, within 

the percentage error of 1% agrees with the Borel summation and our own estimate for v. 
Therefore, our estimate for a may be considered as satisfactory. 

The critical index 8, with the n = —2 limit violated by the Wilson e— expansion was 
estimated in Section 6. We found that S(n = 0) = 4.882. Theoretical field Borel summation 
data are not available for comparison. From our estimate 773 = 0.023 and the scaling relation 
5 = j^jj, we obtain 5 = 4.865, and within the percentage error of 0.349% both our estimates 
agree. 

We believe that both good and reliable estimates can be obtained only from expansions 
possessing correct limits by design. Examples of such behavior are given by the e— expansion 
and field theory expansion in powers of an interaction constant. 

IX. APPLICATION TO e-EXPANSION 

In accordance with all our previous remarks, concerning the importance of the correct 
limits at n = — 2 and n — > 00, we considered the well-known Wilson e— expansion [1] and 
found that it does not obey these limits for the critical indices 5, 7, a and v. 

We have seen in Section 6, that for the critical index 8 both limits are violated with 
a percentage error equal to —10%. The critical index 7 (see (83)) does obey the n = — 2 
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limit, but at n — > oo, 7 = 1.75, and the percentage error equals —12.5%, when compared 
with the exact 7 = 2. The critical index a (see (84), obeys the n = — 2 limit, but at 
n — > 00, a = —0.75, instead of the exact ct = —1, with the percentage error —25%. The 
critical index v (see (76)) obeys the n = — 2 limit, but at n — > 00, z/ = 0.875, with the 
percentage error —14.286% when compared to the exact u = 1. Clearly these discrepancies 
should lead to an uncontrolled error within the physical region n — 0, 1, 2, 3. 

Fortunately, the last two indices, rj and (3 (see (82) and (79)) do obey the n = — 2 and 
n — > 00 limits: 77 = at n — —2, n — > 00 and /3 = | at n = —2, and /? = | at n — > 00. 

Compare, in the physical region, the values of 5, 7, a and z/ obtained from r\ and /5 by 
means of the scaling laws with those obtained by the direct use of the Wilson e— expansion. 
The percentage deviation from the initial Wilson expansion, for the index 5 is (in this 
Section, in order not to cause confusion, we use the word "error" instead of the letter "e" 
in the formula (59)): 

E(5,n = 0) = 7.7%, E(5, n = 1) = 7.21%, 
E(5, n = 2) = 7.06%, E(5, n = 3) = 7.082% 

For the index 7 : 

E( 1 ,n= 0) =0.7%, E{ 1 ,n = I) = 1.97%, 
E(rf, n = 2)= 3.27%, £(7, n = 3) = 4.49%, 
the error increases considerably with n, while at n = the correct limit n = — 2 still favorably 
influences the results. For the index a we have 

E(a, n = 0) = -0.214%, £?(«, n = 1) = -20.36%, 
E(a,n = 2) = 168.95%, E(a,n = 3) = 52.19%, 
again being negligible for n = but growing with n. For the index v we get 

£(z/, n = 0) = 1.67%, E(is, n = 1) = 3.13%, 
£(z/, n = 2) = 4.52%, E(y, n = 3) = 5.75%. 
We see that only at n = the results possess a reasonable accuracy, and the quality of 

the Wilson e— expansion is good enough to reach, e.g. by means of n + 2— expansion (see 

Section 8), the quality of the best known estimates. 
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We conclude that all attempts to improve the critical indices 8, 7, a, v for n = 1,2,3 
directly from the Wilson e— expansion will leave us with an uncontrollable error. It seems 
reasonable to renormalize self-similarly only the critical indices r\ and f3 possessing correct 
limiting values and to calculate all other indices form the scaling laws. 

For the index r\ (see (82)), using e as a parameter for renormalization, the following 
approximations are available 

772(e) = a 2 (n)e 2 , a 2 (n) = ^JIjt, 

%(e) = ?72(e) + a 3 (n)e 3 , a 3 (n) = ( 2 72 + 56n - n 2 ). 

The multiplier 777,3(71, s) at e = 1 is always minimal at s — > 00, therefore 

% *(6)=a 2 (n)6 2 exp{^ e r}. (85) 

For i* = 1 we obtain 

77*(n = 0) = 0.045, r? 3 *(n = 1) = 0.051, 
r}*(n = 2) = 0.052, ?7*(n = 3) = 0.05. 
These values are probably too large. The situation may be improved if (85) is optimized 
using the knowledge of the exact 77 = 0.2083 for e = 2, n = (2-d random walks) [27]. 
Setting 773*(e = 2, n = 0,t*) = 0.2083 we find that t* = 0.567 and 

r7*( n = 0) = 0.029, r]*(n = 1) = 0.033, 
77* (n = 2) = 0.034, 77^ (n = 3) = 0.034, 
in a perfect agreement with the best estimates by the Borel summation [20]. If t* is optimized 
from the knowledge of the 2d Ising model exact rj = 0.25, then similarly, t* = 0.603 and 

r)*( n = 0) = 0.03, r/*(n = 1) = 0.034, 
r)*(n = 2) = 0.035, ^(n = 3) = 0.035. 
There is also another way to get information from the e— expansion for the index 77. 
Proceed similarly to the case of ^-—expansion, and re-expand (82) in powers of y — (at 
e = 1). Then, up to the second order in y, 

3 n 2 
i]^ -y + 9y 
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and r/2 = §yexp(24y), with 

rj*( n = 0) = 0.025, r]*(n = 1) = 0.034, 
r)*( n = 2) = 0.039, r£(n = 3) = 0.042, 
still in a good agreement with a set of data available for the critical index 77 [21]. 
For the critical index f3, Eq.(79) defines the following set of approximations: 

A)( e ) = a o, «o = 1, 
0i(e) = A)( e ) + ai(n)e, ai(ra) = -j^gj, 
/3 2 (e) = p 1 (e) + a 2 (n)e 2 , a 2 (n) = ^gg^ - 

The multiplier 

ai(n) 1 + s a 2 (n) 2 + s 2 

m 2 m, s) = 1 H e H e 

a s a s 

is equal to zero at s — — ggaia an( j the evolution integral gives 

^ a2+ai+ao' ° ° 

= ^ — 2a2 4^i) 1 . • ( 86 ) 

From here, at e = 1, 

/3*(n = 0) = 0.313, (3* 2 (n = 1) = 0.333, 
0*(n = 2) = 0.35, (3*(n = 3) = 0.364. 
The formula (86) is applicable only up to n = 42, where s — > and (86) becomes undefined. 
For n > 42, the exponential summation is optimal: 

%(e)=/3 1 (e)e X p(^e), n>42. 
ao 

For n — 2, 3, f3 2 almost coincides with the Borel summation values 0.346 ± 2 and 0.365 ± 2, 
respectively, being larger for n — 0, 1, where the Borel summation gives 0.302 ± 15, and 
0.325 ± 1 respectively. 

The e— expansion was obtained also from the theoretical field approach up to the fifth 
order in e [28]. For the critical index v, up to the second order in e, one has 

V = v^ 1 = a a (n) + ai(n)e + a 2 (n)e 2 , 

(87) 

a a (n) = 2, ai(n) = -^, a 2 (n) = -^fg^, 
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the limits n = — 2 and n — > oo being satisfied, in distinction from the original e— expansion. 
We use the following approximations to the quantity V = —V + a (n) : 

Vi(e) = -ai(n)e, 

(88) 

^(e) = — ai(n)e - a 2 (n)e 2 . 

All terms in (88) are positive and, the optimal renormalization corresponds to s — > oo. The 
evolution integral can be readily calculated giving 

V*(e) = -ai(n)eexp{^re **}■ (89) 

ai(n) 

We found that at t* = 1, e = 1, 

!/*(n = 0) = 0.607, u*(n = 1) = 0.655, 
z/ 2 *( n = 2) = 0.698, i/*(n = 3) = 0.736. 
Similarly to the case of the index 77, the self-similarly renormalized e— expansion for 
v overestimates the critical index, as may be seen from the comparison with other results 
[21]. Let us optimize the expression for v\ using the knowledge of the 2d Ising v — 1 [29]. 
From the condition v^ie = 2,n = l,t*) = 1 we find t* = 0.576 and the optimized values 

u*( n = 0) = 0.59, u*(n = 1) = 0.628, 
z/ 2 *(n = 2) = 0.662, i/*(n = 3) = 0.691. 

Compared to the best known calculations of v from the Borel summation and similar 
methods, for n — 0, 1 our estimates practically coincide with them, and for n = 2, 3 the 
percentage errors are —0.451% compared to v — 0.665 at n = 2 [30] and 1.003% compared 
tou = 0.698 at n = 3 [30]. 

It is interesting that by a single parameter optimization the index v is improved in the 
whole physical range. It is worth noting that t* used for optimization is about the same 
for i] and v. Evaluation of the critical indices is also possible basing on the field-theoretical 
expansions. Information is obtained in the latter case from the perturbative series directly 
and the results are marginally sensitive to the way in which the position of an infrared stable 
fixed point is determined [20,31]. 
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X. APPLICATION TO THE FIELD THEORY EXPANSION 



Theoretical field approach in the theory of critical phenomena gives, probably, the most 
accurate and consistent estimates for the critical indices r\ and 7 [2]. The analysis of the 
expansions in powers of the interaction constant g (^-expansion) for these indices from the 
viewpoint of the limiting cases n = —2, n — > 00, became possible only when the ^-expansions 
have been written for arbitrary n [30]. By direct inspection of the expressions for 77 and 7 
from [30], we found that n — > 00 limit is obeyed rigorously if g — 1, i.e i] = 0, 7 = 2 and 
the n = — 2 limit is obeyed with very high accuracy for arbitrary g, i.e 77 ~ 7 ~ 1 with the 
error insignificant within the framework of the theory of critical indices. 

The standard approach [20] uses, for computing the infrared stable fixed point g* of the 
beta-function W(g), a complicated Borel summation technique. Then critical indices are 
calculated as 7(<7*), v(9*)- This approach requires a number of terms in the expansions. We 
apply below the self-similar renormalization to only the initial three terms in the expressions 
for W(g), 7 _1 (fi'), i](g) and obtain estimates with an accuracy comparable to the best known 
Borel summation results obtained from all known terms in the expansions. 

We construct the following set of approximations to W(g) using the g— expansion from 
[30]: 

w 2 (g) = -9 + g\ 

W 3 (g) = W 2 {g) - a 3 (n) g\ a 3 (n) = 6^4074o^±p^ 

From the equation W^fiO — f we find the expansion function x(f) = \{1 + y/1 + 4/ ). The 
points of the approximation cascade trajectory are 

y 2 (/) = /, 
y3(f)=y2(f)-fxV)- 

For the velocity field we get 

vs(f) = -^ 3 (/). 
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By analogy with Section 4, the evolution-integral solution for W£ is obtained implicitly from 
the equation 

W df _ 

Jw 2 (g) V 3 (f) 

and the root g* of the equation W^(g, n, t* = 1) = is obtained numerically, as a function 
g* = g*(n). In the physically important cases 

g *(n = 0) = 1.59, 0*(ra = 1) = 1.559, 
g *(n = 2) = 1.524, £*(n = 3) = 1.491. 

At n = —2, g* = 1.599 and at n — > oo, g* = 1. The dependence of g*(n) in the interval 
n = —2, is nonmonotonous, a maximum is reached at n — —1, where g* = 1.61. Our values 
are higher than the results of [30], but remain within the reasonable bounds and show the 
same tendency, at least for n = 0, oo. No data are available for comparison for n = —2, 0. 

For the critical index rj, keeping the starting two terms in powers of g, we can write down 
the following approximations 

rj 2 (g) = b 2 {n)g\ b 2 (n) = 2^2™^ 

%(^)=^(g)+63(^ 3 , b 3 (n) = 0.024^14n» + 0.246^14n + 0.3949440224 > 

and the evolution integral can be readily calculated resulting in 

V*s(g,n) = b 2 {n)g 2 exp{j^-g}, 

and 

v *{n = 0) = 0.027, r/ 3 *(n = 1) = 0.030, 
77* (n = 2) = 0.031, r/|(n = 3) = 0.031. 
These values are practically the same that quoted in [20] where the Borel summation has 
been used: rj(n = 0) = 0.026 ± 3, rj(n = 1) = 0.031 ± 4, rj(n = 2) = 0.033 ± 4, r/(n = 3) = 
0.033 ± 4. At n — 4, r/| decreases to 0.03, the same tendency is seen in the data of [30]. 
For the critical index 7 we keep the three starting terms in powers of g : 

n + 2 n + 2 

G^- 1 = l + ci(n)g-c 2 (n)(7 2 , ci(ra) = ~ 2 ( n + 8 ) > = ~ ( n + 8 ) 2 - ( 92 ) 
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The following approximations to the quantity G = —G + 1 are used: 



The multiplier m 2 (g,n, s) 




s = 0. The evolution integral gives 



G* 2 (g,n) = ci(n)g 



1 



1 + 



C2(n) 
ci(n) 



For 72 we obtain 



7 *(n = 0) = 1.166, 7 2 *(n = 1) = 1.239, 



7 *(n = 2) = 1.305, 72 *(n = 3) = 1.363. 



The percentage errors are e 2 ( 7 (n = 0)) = 0.431%, as compared to 1.161 from the Borel 
summation: zero error as compared to 1.241 ± 2 at n — 1; e 2 ( 7 (n = 2)) = —0.836% 
as compared to 1.316, and e 2 ( 7 (n = 3)) = —1.659%, as compared to 1.386 from the Borel 
summation. In the latter case of n = 3 we also constructed the set of approximations directly 
for the index 7 , expanding (92) in powers of g, up to the third order : 

7 = 1 + d\(n)g + d 2 (n)g 2 , di(n) = -ci(n), d 2 (n) = c\(n) + c 2 (n). (93) 

Approximating 7 = 7 — 1 by 



At t* — 1 we obtain 7 |(n = 3) = 1.37 and the percentage error —1.154%. The same 
procedure applied to n = 0,1,2 always keeps the error e < 1%. These results are quite 
accurate, especially if to remember that only the starting terms were used. Unfortunately, 



7l (n) = d!(n)g, 
= 7i( n ) + d 2 {n)g 2 , 



and finding out that the multiplier m 2 (g,n, s) is minimal at s — > oo, we have 



72 *(n) = di(n)gexp{ 



d 2 (n) 
di(n) 



g t*}. 



(94) 
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the expression (93) does not obey the n — > oo limit with a percentage error —3.7%. An 
optimization of the expression (94) may be carried out requiring that the limit 7 = 2 
at n — > 00 be restored (see also section 5). This aim is achieved with t* = 1.15 and 
72(71 = 3) = 1.375 with a percentage error —0.794%. An effective increase of t* mimics the 
effect of introducing higher terms into consideration. 

XI. CONCLUSION 

We suggested here a new variant in the method of self-similar renormalization permitting 
to find effective sums of asymptotic series. The advantage of the method is that it allows 
to get results exploiting just a few first terms of given series, when no other resummation 
techniques work. At the same time, the accuracy of the results is not worse than that 
reached in other known sophisticated techniques involving about ten terms. In addition, our 
method in the majority of cases, makes it possible to present answers in the form of simple 
analytical expressions that are easy to study for considering their dependencies on various 
parameters, including asymptotic behavior with respect to these parameters. 

The possibility of realizing a renormalization, having in hands only a few terms of a 
series, is due to an algebraic transformation which is equivalent to the effective increase of 
approximation orders. That is why we call this variant the algebraic self-similar renormal- 
ization. 

The general idea of the self-similar renormalization [4-8] is to extract the maximum of 
information from the minimal number of terms. Such a minimax criterion, certainly, can be 
followed only with the help of additional functions making the convergence as fast as possible. 
These functions are called control or governing functions, and they play the same role as the 
control functions in the optimal control theory. In the algebraic self-similar renormalization, 
the control functions are introduced into powers of an algebraic transformation. The choice 
of these control functions is based on the principle of maximal stability by minimizing the 
absolute value of mapping multipliers. 
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We illustrated the effectiveness of our approach by renormalizing divergent series in 
the theory of critical phenomena. Doing this, we specially restricted ourselves from using 
many terms of perturbative series, which are sometimes known-This is to emphasize that 
our approach is effective when, really, only a minimal information is available. Dealing 
with higher-order terms needs a multiple repetition of the renormalization procedure. This 
requires a slight generalization of the technique and much more place for presentation. This 
multiple renormalization will be the subject of a separate publication. 

The method suggested is quite general and can be applied to arbitrary divergent series. 
The choice of examples from the theory of critical phenomena owes to the common interest to 
this problem. We wanted also to stress that even for the problem, where so much is known, 
there are ways of improvement by first, obtaining the results much easier, second, deriving 
analytical formulas, not involving heavy numerical calculations, and, finally, by restoring 
correct asymptotic behavior with respect to physical parameters, such as the number of 
components. 

Another important message which we would like to bring up to readers is that one should 
not be afraid of simple perturbative series that, being divergent, seem, at the first glance, 
to be senseless: Even a seemingly bad perturbative series contains quite a lot of useful 
information which can be extracted by means of an efficient renormalization procedure. We 
hope we were able to convince the reader that the algebraic self-similar renormalization 
suggested can be such a tool for extracting quite accurate information even from bad and 
short series. 

This work was supported by a grant from the National Science and Technology Deve- 
polement Council of Brazil. 
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